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I . INTRODUCTION 




It is now well accepted that the heating which maintains the 
corona at a much higher temperature (> 10 K) than the photosphere 
and chromosphere is due to the dissipation of mechanical waves which 
are generated by the "turbulent" motions in the convection layer just 
below the photosphere (Kuperus, 1969) • It is also fairly well esta- 
blished that these waves are compression waves with periods in the 
range 200 to 300 sec (Uchida, 1967; Kopp, 1968; Moore, 1972 ). However, 
the manner in which these waves dissipate in the corona is not well 
understood. In the first place, even the basic physical mechanism by 
which the waves are damped is not established. Possible damping mecha- 
nisms are shock-wave dissipation (Kopp, 1968), damping by thermal con- 
duction (Sturrock, 1964; Moore, 1972 ), and Landau damping (D’Angelo, 
I 969 ; Moore, 1972). Each of these mechanisms seems to be capable of 
dissipating the compression waves in the corona. Secondly, assuming 
that the waves are dissipated by one of these processes, the theore- 
tically computed rate of heating due to the dissipation at any given 
point in the corona is rather uncertain. For these reasons, it is of 
interest to empirically determine as far as possible the mechanical-wave 
heating in the corona. 

The mechanical-wave dissipation may be characterized by the damping 
length Lp defined by 


1 dF 

T = _ I ™ 

D ~ F dz 9 
m 


(1) 


where is the flux of mechanical energy carried upward by the waves 


and z is the vertical coordinate. The purpose of this paper is to 



obtain an empirical estimate of in the lower corona in quiet 

regions. Our approach is to first empirically determine the relation 
between the dissipation and the temperature structure of the lower 
corona. We then use this relation in conjunction with the observed 
structure to estimate the damping length. 
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II. RELATION OF THE DISSIPATION AND THE TEMPERATURE STRUCTURE 

It is known from XUV- resonance- line observations that between 
the chromosphere and the corona there is a narrow transition region 

4 

in which the temperature increases from about 10 K to the vicinity 

g 

of 10 K in a height range of a few thousand kilometers (Dupree and 
Goldberg, 1967 ) • In the lower corona the temperature increases more 
slowly with height, reaching a flat maximum of about 2 X 10 K at a 
height of several tens of thousands of kilometers (Withbroe, 1970) ♦ 

The temperature level at which the transition region ends and the 
corona begins is somewhat arbitrary, any temperature within about a 
factor of two of 10^ K being reasonable. For the purpose of this 
paper, it is convenient to adopt the 6 X lCr K level as the base of 
the corona . 

We can determine the relation between the dissipation and the 
temperature structure of the lower corona by considering the energy 
balance of the corona as a whole. The heating of the corona is 
balanced by three modes of cooling: downward heat conduction to the 

transition region, radiative cooling and energy loss to the solar wind. 
In the following paragraphs we estimate the empirical rate of energy 
loss from the corona due to each of these modes. 

The temperature gradient in the upper transition region and lower 
corona can be estimated from the observed limb brightening of XUV 
resonance lines emitted from these layers (Withbroe, 1970) . At the 
6 X lCr K level, the empirical temperature gradient obtained by 
Withbroe is 

m , 10 -2- 4 5 ± 0-3 K cm -l . ( 2 ) 

We use the subscript zero (0) to denote the base of the corona (where ■ 
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R 

T = Tq = 6 x lCr K) . At any level in the atmosphere the heat flux 

F is proportional to the temperature gradient at that level: 
c 


F 

c 



( 3 ) 


where H is the thermal conductivity. In the temperature and pressure 
range of the upper transition region and lower corona, the thermal 
conductivity is given by 

,^-5.18 m 2.3 6 -1 , -1 -1 

K = 10 J T J erg sec K cm 

(Moore and Fung, 1971 ) . Combining Equations (2), (3) and (4), we 
obtain 

F = -10^*^ ± 

c o 

for the observed value of the rate of energy loss from the corona by 
downward heat conduction. 

In quiet regions, almost all of the energy radiated from the upper . 
transition region and the corona is contained in XUV emission lines 
(Pottasch, 1964). From identifications of the ions and the transitions 
which produce the lines in the solar XUV spectrum and from estimates 
of the temperatures at which these lines are most likely to be excited, 
Nikolsky ( 1969) has estimated the total radiative energy loss above 

R 

the 6 X lCr K level. For the quiet solar atmosphere, Nikolsky finds 

from the observed spectrum that the radiative cooling rate of the 

5 -2 -1 

corona is less than 10 erg cm sec 

Since the solar wind is highly supersonic at the orbit of the earth, 

the thermal energy density is negligible compared to the kinetic energy 

density of the streaming motion. The energy flux of the solar wind can 

therefore be estimated from the mass density p and velocity v of 

r E E 


-2 -1 

erg cm sec 


(5) 



k 



the solar wind observed at the orbit of the earth. The energy flux 

F at the surface of the sun necessary to maintain the solar wind 
sw 

is given by 

F S„ - e E v E (s v e 2 - (4) ’ (6) 

where GM^/R^j is (very nearly) the gravitational potential energy per 
unit mass at the orbit of the earth, and R„ the radius of the earth's 
orbit. During quiet times, the average number density of ions (mostly 

_ o 

protons) at the earth is typically 10 cm , and the velocity is typi- 
cally kOO km sec * (Ness, 1968) . These values give 

F gw = 8.7 X 10 ^ erg cm 2 sec 1 ( 7 ) 

for the rate of energy loss from the quiet corona to the solar wind. 

The preceding estimates show that the cooling of the corona is 
completely dominated by downward heat conduction to the transition re- 
gion. This requires that almost all of the mechanical energy flux be 
dissipated below the temperature maximum in the lower corona. There- 
fore, (l) the upward mechanical energy flux F which enters the 

m 0 

corona is approximately equal to the downward heat flux F which 

c o 

flows out of the corona: 



( 8 ) 


and (2) at each point in the corona below the temperature maximum, the 
rate of mechanical-wave heating per unit volume is approximately equal 
to the rate of conduction cooling per unit volume: 


dF dF 
m c 

— — — — 

dz dz 


(9) 


From Equations (l), ( 3 ) , (8) and (9), we obtain 
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( 10 ) 


in the corona below the temperature maximum. This is the basic rela- 
tion which connects the wave dissipation with the temperature struc- 
ture of the lower corona. 

At the base of the corona F = F » -F and T = and at 

0 c 0 "o 0 

the temperature maximum, F =0 and T = T . Integration of 

c max 

Equation (10) from the base of the corona to the temperature maximum 
gives 




(ii) 


H is a function of temperature, and may vary with height in the 

corona; but to illustrate the effect of the dissipation on the tem- 
perature of the corona, we temporarily assume for simplicity that both 
the thermal conductivity and the damping length are constant in the 
lower corona. Then Equation (ll) becomes 




(T 


max 



( 12 ) 


Thus, we see (l) that a larger total amount of heating F^ results 
in a higher maximum temperature in the corona as is to be expected, 
but (2) that for a constant total amount of heating F^ , stronger 
damping (shorter damping length) results in a lower maximum tempera- 
ture in the corona. 
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III. ESTIMATE OF THE DAMPING LENGTH 


From Equation (ll) we see that an "average” or representative 
value of the damping length in the lower corona is given by 


J 

J D = t=f V J 

( e o) T t 


max 


H dT . 


(13) 


Using Equation (4) for H, we obtain 


L d = 2.0 X 10 


m3 • 36 / T \3 • 36-1 

-6 max Tn / 0 \ 

T^~\ 1 " It — ) I 

I C Q ) *“ ' max ' J 


cm . 


Since T A = 6 X 10^ K and T ss 2 X 10^ K, (T./T ) 3 ’^ 6 

0 max 3 v 0 max' 

so that we may take 


(14) 


1.8 X 10 2 , 


L d = 2.0 X 10 


T 3.36 

-6 max 


r%r 


cm . 


( 15 ) 


We estimate from Withbroe's (1970) results that the limb brightening 


data determine the value of T ° / 

max 


30 $: 


K) 


with an uncertainty of about 


max 


TO 


1 A 12.9 ± 0.1 „3 -1 2 

10 ^ K erg cm sec . 


( 16 ) 


Withbroe finds that the value of T can be determined from XUV- 

max 

resonance-line intensity ratios, also with an uncertainty of about 30 $: 


T = 10^*3 ± 0,1 K 
max 


(17) 


Using these values in Equation ( 15 ) we obtain 
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(18) 


— Q 

L d = 3.0 X 10 y cm 

with an uncertainty of about b0 % . This estimate of the damping length 
provides an observational criterion which may be used to test theore- 
tical treatments of the mechanical-wave heating in the lower corona. 
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